Charge expulsion and electric field in superconductors 
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The theory of hole superconductivity predicts that when a metal goes superconducting negative 
charge is expelled from its interior towards the surface. As a consequence the superconductor 
in its ground state is predicted to have a non- homogeneous charge distribution and an outward 
pointing electric field in its interior. Here we propose equations to describe the behavior of the 
charge density and electric field in superconductors, and solve them for a spherical geometry. The 
magnitude of the predicted interior electric field depends on superconducting parameters such as 
the condensation energy and the London penetration depth and is found to be of order 10^1^/ cm. A 
physical interpretation of the result is given. It is predicted that for small superconducting bodies 
(compared to the penetration depth) an electric field outside the superconductor should result from 
this physics. This may explain a recent experimental observation in Nb metal clusters. 

PACS numbers: 



I. INTRODUCTION 

In the currently accepted understanding of supercon- 
ductivity no electric fields exist in superconductors in 
the absence of electric current (a 'Bernoulh potential' 
is expected to exist in the presence of a non-uniform 
supercurrentQ). The London brothers originally pro- 
posed a set of equations to describe superconductivity 
that allowed for the presence of an electric field within 
a penetration depth of the surface of a superconductor. 
However H. London failed to detect such an electric field 
experimentally 0. Having confidence in his brother's ex- 
perimental result, F. London modified the theory and 
discarded one of their original equations, and as a conse- 
quence the possibility of an electric field in superconduc- 
tors is no longer discussed in London's definitive workQ. 

As is well known, normal metals allow for the existence 
of magnetic fields but no electric fields in their interior. 
Superconductors (type I) do not allow magnetic fields in 
their interior, but magnetic fields can exist within a pen- 
etration depth of the surface. Based on the theory of 
hole superconductivity!^ ^ we have recently proposed 
that when a metal goes superconducting negative charge 
is expelled from its interior towards the surface!^. Here 
we discuss the equations governing the charge and elec- 
tric field distribution in superconductors resulting from 
this physics. The possibility that a superconductor may 
have an electric field in its interior has not been dis- 
cussed in other theoretical frameworks to our knowledge; 
it is however a necessary consequence of the fundamental 
electron- hole asymmetry of condensed matter and the 
resulting 'giant atom' description of superconductors 
that results from the theory of hole superconductivity. 



II. ENERGETICS 




FIG. 1: Schematic picture of a spherical superconducting 
body. Negative charge is expelled from the bulk to the surface. 
As a consequence, an outward-pointing electric field exists in 
the interior. 



ant atom, with a higher density of negative charge near 
the surface and higher density of positive charge in the 
interior Consider a superconducting sphere of radius 
R much larger than its London penetration depth A^. 
Assume an amount of charge q moves from the interior 
to the surface, resulting in a uniform charge distribution 
Po in the interior: 
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As we discuss later, q is negative and resides within a 
London penetration depth of the surface, with volume 
charge density 



The qualitative picture of a superconductor proposed 
in refs. 0, is shown in figure 1: it looks like a gi- 
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For Al < < i? we can also think of this charge as a 'surface 
charge density' 
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This 'surface charge density' is not confined within a 
Thomas Fermi screening length of the surface as excess 
charge would be in ordinary metals, but within the much 
thicker surface layer defined by the penetration depth. 
The electric field in the interior is given by 
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and the inhomogeneous charge distribution gives rise to 
a cost in Coulomb energy 
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under the assumption Xl « R- For a spherical geom- 
etry no electric field exists outside the superconductor 
(assuming charge neutrality). 

In the superconducting state, the energy per unit vol- 
ume is lowered by the superconducting condensation en- 
ergy. The theory of hole superconductivity predicts that 
the charge expulsion describe above is a necessary conse- 
quence of superconductivity. Hence the Coulomb energy 
cost Eq. (5) is balanced by an energy gain due to the 
charge expulsion, e per unit volume, which is related but 
not necessarily equal to the condensation energy of the 
superconductor. Equating this energy gain to the elec- 
trostatic energy cost, 
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yields for the 'surface charge density' 
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independent of the sphere's radius. This indicates that if 
the physics discussed here indeed exists in superconduc- 
tors, the charge density in the surface layer p- is an in- 
tensive quantity determined by intrinsic properties of the 
superconductor, in particular the penetration depth and 
the thermodynamic critical field He (which gives the con- 
densation energy). Instead, the interior positive charge 
density po is not intensive but decreases with increasing 
radius R. The electric field close to the surface, given 
by Eq. (4) for r ~ i?, is also independent of the body's 
dimensions. 

From Eqs. (3) and (7) the volume charge density in 
the surface layer is given by 



with m,. the free electron mass, e the electron charge 
(negative) and the superfluid density per unit vol- 
ume, and replacing e in terms of the energy gained per 
superfluid electron e yields 



e = 



and Eq. (8) becomes 
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relating the negative excess charge density in the surface 
layer to the superfluid density Ug and the ratio of the 
condensation energy of an electron e to its rest energy, 
showing that the excess charge density is a small fraction 
of the superfluid density. The 'surface charge density' 
a is given by Eq. (3), and the maximum electric field 
attained within a penetration length of the surface of the 
superconductor is 
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and is also independent of the sample's dimensions. 

As an example appropriate for high cuprates we take 
Xl = 2000 A, yielding from Eq. (9) a superfluid density 
rig = 7.1 X 10~^ electrons /A^. Assuming a condensation 
energy of lOOpeV per unit cell and a unit cell volume 
of 40A^ yields e = 3.5meV per superfluid electron, and 
P- = —1.5 X 10~*nse so there is roughly 1 extra electron 
per 10, 000 superfluid electrons in the outer layer. The 
'surface charge density' is cr = —3.5 x lO^'^C/cm^, and 
the maximum electric field near the surface is E^^x = 
3.9 X \(fV/cm. The potential difference between the cen- 
ter of the sphere and a point within a penetration depth 
of the surface is A]/ = _E„iarE-R/2, so approximately 4 
million Volt for a sample of 1cm radius. 

As another example, for Nb the thermodynamic crit- 
ical field is He — 1980G, so the condensation energy is 
e = 1.56x lO^ergs/cm"^. From Eq. (7) the 'surface charge 
density' is cr = 6.8 x 10~^C/cm^ and taking Xl = 400A 
for the penetration depth yields p_ = 0.017C/cto^ for the 
surface layer volume charge density, corresponding to ap- 
proximately 2 excess electrons per million Nh atoms in 
the outer layer. The maximum electric field near the sur- 
face from Eq. (12) is E^ax = 0.77x lO^F/cm. This is not 
very different from the estimate for the high Tg case be- 
cause the efii'ects of a smaller condensation energy and a 
smaller London penetration depth partially compensate 
each other. 
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Using the expression for the London penetration depth 
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III. ELECTRIC FIELD EQUATIONS 

We start with the London equation for the supercur- 
rent 
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with A the magnetic vector potential. Following the 
Londons0 we assume that A satisfies the Lorenz gauge 
condition 
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with <j) the electric potential. The electric field is given 

by 



E = -V0- 
Using the continuity equation 
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and applying the divergence operator to both sides of Eq. 
(13) and using the gauge condition Eq. (14) yields 
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hence 
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The London brothers postulated Eq. (18) with (l)o{'r) = 
as a possible equation applicable to superconductors 0- 
Under stationary conditions, 



(t){f) = -47rAip(f) 
and using Maxwell's equation 

V • E{f} = Anpir) 
Eqs. (15) and (19) yield 

p(r) = Aiv2p(rO 
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which imply that a non-zero charge density and an elec- 
tric field can exist within a penetration depth of the sur- 
face of the superconductor but not in the interior. H. 
London attempted to measure the electric field predicted 
by Eq. (21) by measuring the change in capacitance of a 
capacitor with superconducting plates when it is cooled 
below Ted. However he failed to detect any effect. 

Here we propose that Equation (18) is valid for super- 
conductors with a non-zero (poir) resulting from a uni- 
form charge density po deep in the interior of the su- 
perconductor, as discussed in the previous section. Eqs. 
(21) and (18) then become 



p(r) = Po + AivV(r) 
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with Ea[r) = — V0o(r) and V ■ Eo[r) = 47rpo- We pro- 
pose that these equations hold for the interior of super- 
conductors of arbitrary shape in the absence of magnetic 
fields and electric currents, with po a- positive number 
determined by the microscopic parameters of the super- 
conductor as well as the geometry, as discussed in the 
previous section. The potential obeys 



0o(f) + Aiv20(rO 



in the interior, and 
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in the exterior of the body. Furthermore we assume no 
surface charges can exist in the superconductor, hence 
that both (j) and its normal derivative d(j)/dn are contin- 
uous on the surface of the body, as are of course and 
d(j)Q/dn . 

Note that the existence of an electric field in the in- 
terior of the superconductor does neither imply a time 
variation of the current, as one would expect in a 'per- 
fect conductor', nor even the existence of a current as 
in an ordinary metal. Taking the time derivative of Eq. 
(13) and using the Maxwell equation (15) yields 
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so that in a stationary situation there can be an electric 
field that is derivable from a potential {E — — V0) and it 
does not lead to a time- varying supercurrent. Whether 
a stationary supercurrent exists or not depends on the 
magnetic vector potential A through Eq. (13) and not 
on the electric field. 

In a spherical geometry, the solution of Eq. (22a) is 
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and for charge neutrality the constant k is given by 
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FIG. 2: Radial dependence of charge density (a) and electric 
field (b) for RJXl = 0.5,1,3.5,10 and 100 (some numbers 
next to the lines). The electric field increases monotonically 
with increasing R/Xl- The charge density at the origin is 
smallest for R/Xl ~ 100, and decreases with decreasing R/Xl 
for R/Xl < 3.5. 



/(x) = a; cosha; — sinhx (26b) 
For i? >> Al the charge density at the surface is 



piR) = Po(i 
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in accordance with Eqs. (1), (2), and for r << R — Al, 
p{r) ^ Po- The electric field is given by 
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It goes to zero at r = i? and is maximum at a distance 
of a penetration depth inside the surface, given by 
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in agreement with Eq. (12) for p{R) = p-, a.s expected. 

Figure 2 shows the radial dependence of the charge 
density and electric field for various values of R/Xl- 
When the dimensions of the sample are much larger than 
the penetration depth there is a small uniform positive 
charge density in the interior and the electric field in- 
creases linearly towards the surface. The compensating 



negative charge resides within a penetration length of 
the surface and the electric field drops to zero in that 
region. As the dimensions of the sample decrease the 
positive charge density in the interior increases and is 
maximum for R/Xl ^ 3.5, where po ~ — 0.56p_, and 
decreases again for even smaller samples. This is under 
the assumption that superconducting properties are not 
affected by the size of the sample, i.e. p_ and Al are 
assumed independent of R. In the limit R/Xl << 1 the 
charge density is given by 
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so that it changes sign at r/R — 0.77, and the electric 
field by 



(31) 



IV. MICROSCOPIC INTERPRETATION 

The existence of an electric field in the interior of a 
superconductor may seem surprising. In normal metals, 
electric fields cannot exist in the interior because the mo- 
bile electronic charge will rearrange to nullify such field. 
How can the existence of an electric field in a supercon- 
ductor be understood? 

Consider Eq. (24). It states that no electric current 
will be generated in the presence of an interior electric 
field, provided the field derives from a potential (j). Now 
an electric field exerts a force on electrons, why is it that 
electrons in the superconductor do not get accelerated by 
the interior electric field and generate an electric current? 

The Lorentz force equation in the absence of a mag- 
netic field 
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indicates that the velocity of an electron will change in 
the presence of an electric field. From the relation be- 
tween velocity and current 



J — en^v 



(33) 



one may be tempted to conclude that Eqs. (32) and (33) 
are incompatible with Eq. (24). However this is not the 
case, because Eq. (24) is a local time derivative and Eq. 
(32) is a total time derivative. Using the relationflTj 

dv dv -tv'^ ^ „ . ^, 

-.- + V--.x(Vx.) (34) 

we conclude that a stationary situation with an electric 
field in the interior of the superconductor is indeed pos- 
sible, provided that 
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If we assume that the superfluid electrons traverse cir- 
cular orbits with speed v{r) we have 



\7— w X (V X -y) 
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and the speed is obtained from Eqs. (35) and (28) as 
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which of course simply describes the balance between 
centripetal acceleration and electric Lorentz force 
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as appropriate for a circular orbit. Deep in the interior 
of the superconductor Eq. (37) describes a 'rigid rota- 
tion' of the superfluid, v(r) = cur, with constant angular 
velocity 
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As we approach the surface the angular velocity decreases 
to zero due to the screening of the electric field by the 
negative charge in the surface layer. 

We conclude from this interpretation that the 'rigidity' 
of the wave function 4] that allows for an electric field in 
the interior arises from the fact that superfluid electrons 
traverse macroscopic orbits at high speed; a radial elec- 
tric field will not change the radial speed of electrons for 
the same reason that an orbiting planet is not attracted 
towards the sun: the radial force provides the centripetal 
acceleration for the orbital motion. 

However one may also ask: if there is a finite velocity 
as given by Eq. (37), shouldn't there be a finite current 
and an associated magnetic field? In fact Eq. (33) is 
an oversimplification, which assumes that all electrons at 
position r move with the same velocity v. However elec- 
trons also have a spin quantum number. As the simplest 
resolution to this question we may assume that the ve- 
locity field for spin direction a (with respect to a chosen 
quantization axis) obeys 



Va{r) = -V-cr{r) 



(40) 



with Va- given by Eq. (37). In that case the total charge 
current 



J = e^(^/| +vi) 



(41) 



will be zero and no magnetic field will be generated, how- 
ever there will be a spin current 



Jspin = eY(wt - vi) 



(42) 



which is not zero. The direction of the spin current would 
be determined by spin-orbit coupling as discussed in Ref. 



A different simple solution of Eq. (35) is obtained un- 
der the assumption V x w = 0. In that case, we have 
simply 
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and the velocity profile is readily obtained by integration 
of Eq. (28) as 
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with k an integration constant. Deep in the interior of 
the superconductor we have (for R/Xl >> 1) 
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Eq. (44b) has a simple physical interpretation. The po- 
tential in which the electron moves in the interior of the 
superconductor is a harmonic oscillator potential: 



U{r) 



(45) 



where (f)Q{r) follows from Eq. (23b). The total energy of 
the electron is sum of kinetic and potential energy: 



Et. 



(46) 



and is a constant of motion. Indeed, Eq. (46) is 
the same as Eq. (44b), with the integration constant 
k = 2Etotlme- The allowed values of k are such that 
the maximum elongation r^ax for which v goes to zero 
obeys rmax < R- When k is such that Vmax approaches 
R, the motion is slightly faster than in the harmonic os- 
cillator potential near the region of maximum elongation 
(as given by Eq. (44a)), due to screening by the nega- 
tive charge near the surface. For small R/\l instead the 
velocity is very different than Eq. (44b) and is given by 
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Note also that an electron that moves in a harmonic 
oscillator potential does not give rise to uniform charge 
density, but rather the electron density and hence nega- 
tive charge is largest near the region of maximum clas- 
sical elongation. This is consistent with the assumption 
that the charge density in the superconductor is positive 
in the interior and negative near the surface. Similarly, 
circular 'Bohr orbits' in a harmonic oscillator potential 
'bunch up' for large radius as discussed in Ref. Q giving 
rise to larger negative charge density near the surface. 

The assumption V x w = implies linear motion where 
the electrons oscillate through the origin. Relaxing this 
condition many other trajectories become possible. Deep 
in the interior the motion is simply described by inde- 
pendent harmonic motion in the three directions with 
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frequency given by Eq. (39), which in general describes 
elhptical orbits. The hnear oscillator and circular motion 
described above are two particular examples of these or- 
bits where it is simple to obtain the corrections due to 
the outer layer of negative charge, as given by Eqs. (44a) 
and (37). 



SUPERCONDUCTORS OF GENERAL 
SHAPE 



From Eq. (22a) we have 

~p{r) = ALV2p(f) 



(48) 



for p{r) = p{r) — pQ. Laue has shown rigurouslv|l2j that 
for a body of arbitrary shape the solution of Eq. (48) 
decays to zero rapidly towards the interior, on the length 
scale of the penetration depth A^. Hence we can con- 
clude that quite generally for superconductors of dimen- 
sions much larger than the charge distribution in the 
interior is constant and equal to pQ, and that the expelled 
negative charge resides within A^, of the surface. 
Similarly using that 



V X (V X £;o) = V(V • Eo) 



(49) 



as well as V X .Eo = and \7 ■ Eq = constant we conclude 
that \/^Eo = 0, hence from Eq. (22b) 



V^{E- 



Eo)^^iE 



En 



(50) 



This equation is the same as the London equation satis- 
fied by the magnetic field in non-rotating superconduc- 
tors, giving rise to the Meissner effect. Again, Laue 
showed [l^ that the vector quantity satisfing Eq. (50) 
decays rapidly to zero over a distance of a penetration 
depth from the surface. Hence we conclude that for su- 
perconductors of arbitrary shape away form the surface 
the electric field is given by Eq originating in the uniform 
charge distribution po- For example, for an ellipsoid of 
revolution with symmetry axis along the z direction 



z 



= 1 



the electric field is given bv[T^ 

(Eol 



gVrpoO-p)^ 



4 X 



26 



(51) 



(52a) 



(52b) 



with p = b/a. Eq. (52) is approximately valid for 0.8 < 
p < 5. In particular, the surface of a uniformly charged 
ellipsoid is not an equipotential surface, and this is of 
course also the case for a body of arbitrary non-spherical 
shape. 



In the presence of the expelled negative charge we ar- 
gue that in general the surface of the body will still not 
be an equipotential surface, and hence that electric fields 
will exist outside the superconductor. This is seen as 
follows: we define 



0(r) = 0(f) - 00 (r) 4- AnXlpo 
which obeys the differential equation 
0(f) = Aiv20(f) 
inside the body, and 

y^d,(r) = 



(53) 



(54a) 



(54b) 



outside. Furthermore, both and its normal deriva- 
tive d(f>/dn are continuous on the surface of the body. 
The function satisfying the differential equation (54) 
is uniquely determined by specifying either <f> or d(j)/dn 
on the surface of the body. For Eq. (54b) this is well 
known 13], for Eq. (54a) it follows similarly from Green's 
theorem since for any two solutions of Eq. (54a) with the 
same boundary values their difference U{r) satisfies 



d?r[\VU\ 



El- 



s on 



(55) 



hence U is identically zero in the interior. If no electric 
field existed outside we would have = constant and 
dcj)/dn = on the surface of the body, and the function 
would be overdetermined, i.e. there is in general no 
solution to the differential equation (54) with both and 
d<j)/dn specified on a closed surface. 

To solve the problem in general numerically we suggest 
the following procedure. Assume an initial guess for 0(f) 
on the surface (for example constant). Solve the Dirichlet 
problem for inside and outside, and find 90/9n at the 
surface coming from the inside and the outside, which in 
general will not match. Hence d4>/dn is discontinuous 
at the surface, which implies the presence of a surface 
charge density 

-J_r^^ 1 (^f.\ 

^ . Y ^ jinside ) outside] \^^) 

47r on on 

Since we assume no such surface charge is possible in su- 
perconductors it indicates that the initial guess for 0(f) 
on the surface was incorrect. We next compute the aver- 
age of 90/9n inside and outside the surface, and solve the 
corresponding Neumann problem for that boundary con- 
dition inside and outside. The resulting 0(f) on the sur- 
face will in general be discontinuous, so we calculate its 
average coming from inside and outside and solve again 
the Dirichlet problem with the new boundary condition. 
Thus solving a sequence of alternating Dirichlet and Neu- 
mann problems the solution should converge to a unique 
solution with both and d(j)/dn continuous at the sur- 
face. 
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These considerations also indicate that for supercon- 
ductors of dimensions much larger than the penetration 
depth the negative charge near the surface will arrange 
laterally so as to nullify any tangential electric field and 
hence that the surface of the body will be very nearly an 
equipotential surface. This is because the surface charge 
density Eq. (56) required for it can be achieved by ar- 
ranging the charge in the surface layer of thickness Xl 
appropriately. In other words, the solution to the differ- 
ential equations where the surface charge a in Eq. (56) 
is spread out over the layer of thickness and </> is very 
nearly constant on the surface will be the unique solu- 
tion for a body of dimensions much larger than Al . If </) 
is constant on the surface and the body is charge neutral, 
no electric field lines can exist outside: for any exterior 
electric field line starting on the surface would have to 
return to another point on the surface, and the integral 
of the electric field along that line would yield a differ- 
ence in potential at the initial and final points, both on 
the surface. Hence we conclude that for a charge neutral 
superconductor of large dimensions essentially no electric 
field lines will exist outside the superconductor, and it is 
not possible to detect the non-uniform charge distribu- 
tion in the interior by measurements outside the super- 
conductor. 

Instead, for superconductors of dimensions compara- 
ble to or smaller than the penetration depth the nega- 
tive charge distribution cannot be regarded as a 'surface 
charge', and it will not cancel the potential differences 
at different points of the surface originating in the posi- 
tive charge. Hence the electric field will 'leak out' to the 
exterior of the superconductor and become observable. 

The simplest shape that can give rise to such effect is 
an ellipsoid of revolution, described by Eq. (51). For 
a uniformly charged ellipsoid with charge density po the 
dipole moment (with respect to the center of mass) is 
zero, but the quadrupole moment is given by 

Q = j d^rpir){3z^ - r^) = ^TvaH{b'' - a^)po (57) 

For the charge neutral superconductor, if negative charge 
is expelled towards the surface the charge distribution 
looks qualitatively as in Figure 3. If the dimensions of 
the body are not much larger than the penetration depth 
it is plausible to assume that the negative charge will not 
rearrange laterally but instead will remain approximately 
uniform. It is clear from Fig. 3 that the negative charge 
will give a larger contribution to the quadrupole moment 
than the positive charge. Consequently this will give rise 
to a net negative quadrupole moment for a prolate ellip- 
soid {b > a with b in the z direction) , and a net positive 
quadrupole moment for an oblate ellipsoid (b < a). 

For body dimensions comparable but somewhat larger 
than the penetration depth we can estimate the 
quadrupole moment assuming charge density within 
distance Al of the surface and po in the interior. The net 
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Oblate 




Q>0 



FIG. 3: Charge density in ellipsoids of revolution of dimen- 
sions comparable to the penetration depth (schematic) . There 
will be a net quadrupole moment (negative for prolate, pos- 
itive for oblate) giving rise to an electric field outside the 

superconductor whose magnitude decays as the fourth power 
of the distance to the center of the ellipsoid. 



quadrupole moment that results is 
16 



Q = ^7ra^6(6 - a)\LP- 



(58) 



so that it is negative (positive) for a prolate (oblate) el- 
lipsoid. The electric field at distance r from the center 
outside the body is of order E ~ Q/r^ and should be 
measurable. 

For a body of dimensions much smaller than the 
London penetration depth we can estimate the net 
quadrupole moment as follows. The charge density Eq. 
(30) becomes for the ellipsoid 



(59) 



to ensure charge neutrality, and calculating the integral 
Eq. (57) over the volume of the ellipsoid yields 



Q 



1575 



TTa^b{b - a) (3a^ + 'iaH + Aab"^ + 46^) (60) 



For small distortion from sphericity in terms of i? ~ a ~ 6 



Q = -^7ra^6(6 - o)Al x 



2 PqR' 
15 Ai 



(61) 



which is similar to Eq. (58) for ~ Aj,- If we relate po 
to p_ using Eq. (59) with p_ = p{r = R) Eqs. (61) and 
(58) differ by a factor 4/3. 
For body shapes with less symmetry than the ellipsoid, 

a dipole moment also becomes possible. As an example, 
consider the body of revolution shown in Figure 4, with 
radius 



r{z) = a + {b — a) 



h 



(62) 
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FIG. 4: Charge density in a body of low enough symmetry 
that it can give rise to an electric dipole moment (schematic). 
For dimensions of the body comparable to or smaller than the 
penetration depth there will be a net electric dipole moment 
along the vertical direction giving rise to a dipolar electric 
field outside the superconductor. 



and < z < h. For a uniformly charged such body 
with charge density p the dipole moment relative to (r = 
0, z = h/2) is given by 



P 



12 ^ 



(63) 



We assume again that the expelled negative charge does 
not rearrange laterally for body dimensions not much 
larger than the penetration depth. The net dipole mo- 
ment is then given by 



p^-p^h[h{b'-a')-h'{b 



i"^ + ab + b'^ 
'^ + a'b' + b'^^ 



(64) 



with h' = h- 2\l, a' = t^Xl) - \l, b' = r{h - \l) - Xl, 
and points along the z direction. The sign of p is not 
simply determined, e.g. for < a it is generally positive 
but can become negative for large h. 



VI. DISCUSSION 

We have proposed that the ground state charge distri- 
bution in superconductors is different from what is con- 
ventionally assumed. Rather than being uniform and lo- 
cally charge neutral as normal metals, we propose that 
the charge distribution in superconductors is inhomoge- 
neous, with more negative charge near the surface and 
more positive charge in the interior, reflecting the funda- 
mental charge asymmetry of matter. 

The charge inhomogeneity is a consequence of the fact 
that the superfluid electrons are highly mobile and have 
quantum-mechanical coherence over the macroscopic di- 
mensions of the sample. These electrons have 'un- 
dressed' from the electron-ion and the electron-electron 



interaction!^ ^| and become completely free-electron- 
like, except for the pairing correlations that bind k f and 
—k i electrons. The quantum-mechanical lowering of ki- 
netic energy of the light electrons due to delocalization 
causes the electronic charge density to be larger near the 
boundaries of the sample, just as electrons do not remain 
confined within the dimensions of the positive nucleus in 
an ordinary atom. 

We argue that the electric field in the interior of super- 
conductors is not screened as in ordinary metals because 
the ground state wave function has 'rigidity', and the 
macroscopic quantum-mechanical coherence prevents lo- 
cal deformations that would screen the electric field. As 
discussed in Sect. 4, an electric field can exist in a super- 
conductor if the superfluid electrons are delocalized over 
macroscopic distances and their velocity has the proper 
gradient. At finite temperatures however there are also 
excited quasiparticles that are not macroscopically co- 
herent, and one may wonder why they do not screen the 
electric field. We suggest that they are unable to screen 
the interior field because they are positively charged [T^ 
and as a consequence they are also pushed out towards 
the surface. 

In the Meissner effect, the supercurrents can shield the 
magnetic field so that it is zero in the interior of the 
superconductor only if the body's dimensions are larger 
than the penetration depth; otherwise, the magnetic field 
penetrates the superconductor. Similarly, we have ar- 
gued that the negative charge can shield the electric field 
so that it is zero in the exterior of the superconductor 
if the body's dimensions are larger than the penetration 
depth; otherwise, the electric field can leak out from the 
superconductor to the exterior, unless the body has per- 
fect spherical symmetry. 

The latter fact should allow for experimental detec- 
tion of these electric fields around small superconduct- 
ing particles of non-spherical shape. Remarkably, very 
recently Moro et al reported detection of spontaneous 
electric dipole moments in Nb clusters 17] at low tem- 
peratures. Their observations that the effect is strongest 
for even number of electrons in the cluster and that it 
occurs only at low temperatures suggests that the effect 
is related to superconductivity as proposed by the 
authors. Moro et al concluded from their observations 
the existence of an internal electric field in the interior 
of the metal cluster of the order of lO^V/ cm, and stress 
that such an internal electric field cannot occur in an or- 
dinary metal and suggests a 'rigidity' of the electronic 
wavefunctions as well as a collective effect. The theory 
discussed here predicts that dipole moments should arise 
from clusters of irregular shape such as shown in Fig. 
4. Experimentally, relation between the measured dipole 
moment and the shape of the cluster has not been exam- 
ined, and it would be interesting to do so to compare with 
predictions of this theory. Concerning the magnitude of 
the expected effect, we estimated in Sect. (1) a maxi- 
mum electric field for Nb comparable to the one inferred 
by Moro et al under the assumption that the energy e in 
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Eq. (7) is the condensation energy. However the magni- 
tude of the electric field will be much smaller for cluster 
dimensions smaller than the penetration depth according 
to the results in Sect. Ill, so that our prediction may be 
inconsistent with the observation unless e is assumed to 
be much larger than the condensation energy. 

For small clusters of more regular shape such as ellip- 
soids of revolution, no dipole moment but a quadrupole 
moment should be observed, as discussed earlier. Exper- 
imental determination of quadrupole moments in super- 
conducting metal clusters has not been reported, and it 
would be interesting to search for this effect to compare 
with the theory discussed here. For example, it has been 
reported that tin clusters adopt prolate geometries 0. If 



so, according to the theory discussed here tin clusters at 
low temperatures should exhibit a negative quadrupole 
moment. 

Further consequences of this physics will be discussed 
in future work. 
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